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Abstract 
The reduction of the impact of carbon dioxide on the environment is a key driver of the energy transition, which will 
progressively reduce the use of fossil fuels and will promote the development of viable technologies. This process will be 
based on renewable energy sources to reduce, to any possible extent, the production and the release of carbon dioxide 
into the atmosphere. CCS (Carbon Capture and Storage) projects are more and more an ineludible stage of this process. 
In the last decade, a significant segment of chemical engineering research has been devoted to the management of 
carbon dioxide release scenarios, which can result both from an accidental loss of containment or from an emergency 
venting of carbon dioxide streams. This gas is asphyxiating and toxic, depending on concentration levels. Consequently, 
acquiring familiarity with its dispersion modelling is a key task for process and process safety engineers. On the other 
hand, carbon dioxide presents a very peculiar behaviour, as it is a dense gas and, below the triple point, its release can 
produce solid formation. The article presents a relatively simple multi-stage validated model, covering the release 
scenario from the source, as a heavy gas, up to the neutral Gaussian dispersion, considering also the very low 
temperatures possibly resulting from the JT (Joule Thomson) effect. Two case studies, carried out through the Megaris 
Platform, illustrate the application of the methodology. It can be useful in the design phase and in the QRA studies, to 
properly size and locate the venting stack and to minimize any health upsets for the operators 
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Introduction 


Carbon Capture and Storage (CCS) is a relatively new process of chemical engineering. Its adoption is dictated by the necessity to 
create a barrier to the migration of carbon dioxide to the atmosphere. Consequently, process segments consisting of significantly high 
concentrations of this gas have been introduced, with the result that accidental or emergency releases may cause the formation of 
toxic clouds at ground level. In this respect, the knowledge of the behaviour of the gas is of paramount importance. The adoption of 
validated software is generally effective in providing the outcome of the release scenarios, but some specific evidence is required to 
carry out a more comprehensive analysis, to completely define the risks mapping. Carbon dioxide behaves as a dense gas, both 
because of its molecular weight, which is much greater than that of air, and because frequently it is released at a temperature much 
lower than atmospheric temperature. Figure 1 shows the release of a mixture of carbon monoxide and carbon dioxide, through a hole, 
showing that CO is dispersed as neutral gas (density as air), whereas carbon dioxide drops down. This is the typical scenario of a 
syngas. 
Single models provide the result of single events, whereas in this case, the following multi-stage cloud evolution is expected: 

© momentum-dominated release from a stack or from a hole 

o  slump-down of the cloud to the ground level 

o dense gas horizontal migration until the bulk-density (air+COz mixture) approaches the density of air 

o neutral dispersion 
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Moving from the release to the passive dispersion, the concentration 
behaviour following the gas impingement with the ground is not described 
by any validated model, as the well-known Britter McQuaid method 
describes the dense gas flow just below a certain concentration threshold. 
The available software, unless it is specifically calibrated, does not cover 
the whole scenario, with the result that the final outcome of the dispersion 


End of momentum jet regime, 
can be uncertain. However, the process safety engineer must make sure \ / 
that the toxic effects are fully identified and properly managed. This article \ th. CO 
shows how can this scenario be modelled using simple validated equations, 6 ỌỌ é 
and, above all, it describes the behaviour of the cloud dispersion 
throughout its entire evolution, introducing some novelty in the 


methodology. Two different case studies illustrate the application of the Fig. 1 Dispersion of CO and C02. 
multi-stage approach. 


1.1 Harmful effects of carbon dioxide 


The toxicity of carbon dioxide is currently quantified through the SLOT and SLOD parameters provided by the British HSE. Notably, 
the probit functions for this gas can be determined based on the definition of the two parameters. 

e SLOT: probability of death for 1% of people 

e SLOD: probability of death 50% of people 


Furthermore, the form of the probit function is: 
Pr=a+bln(C™t) (1) 
where: 


Pr is the Probit function, which is Table 1 - SLOT and SLOD for CO». 
related to probability, C is the substance 


Substance n SLOT (ppm"xmin) SLOD (ppm®xmin) 
Carbon Dioxide 8 1.50x1049 150x104! 


concentration, expressed as ppm, f is the 
exposure time, minutes, a and b are 


unitless coefficients, and n is an exponent defined for the specific substance. HSE data for carbon dioxide are shown in Table 1. 
a and b for any substances can be calculated by writing down the probit function values corresponding to 1% and 50% probabilities 
of death from the following table. 


Table 2 Probit functions versus probabilities correlations. 


% 0 1 2 3 4 5 6 i 8 9 
0 - 2.67 2.95 3.12 3.25 3.36 3.45 3.52 3.59 3.66 
10 3.72 3.77 3.82 3.87 3.92 3.96 4.01 4.05 4.08 4.12 
20 4.16 4.19 4.23 4.26 4.29 4.33 4.36 4.39 4.42 4.12 
30 4.48 4.5 4.53 4.56 4.59 4.61 4.64 4.67 4.69 4.72 
40 4.75 4.77 4.80 4.82 4.85 4.87 4.90 4.92 4.95 4.97 
50 5.00 5.03 5.05 5.08 5.10 5.13 5.15 5.18 5.20 5.50 
60 5.25 5.28 5.31 5:33 5.36 5.39 5.41 5.44 5.47 5.50 
70 5.52 5:55 5.58 5.61 5.64 5.67 5.71 5.74 5.77 5.81 
80 5.84 5.88 5.92 5.95 5.99 6.04 6.08 6.13 6.17 6.23 
90 6.28 6.34 6.41 6.48 6.55 6.64 6.75 6.88 7.05 7.33 

0.00 0.10 0.20 0.30 0.40 0.50 0.60 0.70 0.80 0.90 
99 7.33 T3T. TAI 7.46 7.51 7.58 7.65 7.75 7.88 8.09 

Specifically: 
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2.67=a+bln(C"t) (2) 
5.00=a+bln(C"t) (3) 


Being C” ¢ the SLOT and the SLOD values (toxic loads), respectively. Solving the equation system, it results: 


Table 3 Probit function coefficients for CO2. 


Substance n a b SLOT (ppm®xmin) SLOD (ppm®xmin) 
Carbon Dioxide 8 -90.94 1.01 1.5x1040 1.50x1041 


So, the Probit function for COs; is: 
Pr=-90.4+1.011n(C t) (4) 


Equation 4, through Table 2, will provide the probability of death for CO2 inhalation. Other important parameters are the TLV-TWA 
and the STEL, which are provided by HSE (2020) and are 5000 ppm and 15000 ppm, respectively. 


1 Materials and Methods 
1.1 Models 


The dispersion scenario for a vertically released stream of carbon dioxide has been depicted in Figure 2. The gas leaves a stack of 
height hs, with a mass flow rate Wco2 at temperature T,. A momentum-driven initial rise 4h is experienced, after which the flow 
slumps to the ground, due to its molecular weight that is significantly greater than that of air, and, eventually, due to its temperature, 
which, following a depressurisation expansion, might result in lower than ambient temperature. 


air borne plume slumping plume 


plume 
rise 


plume 


touch-down 
Tamb» Weo2 +Wair(Xem) 


> + XL-BMQ > <4 > 


y * Xid XBMma 
HMP Regime Laminar Regime BMQ Regime 


Fig. 2 Dispersion Scenario. 


Four sequential fluid dynamic stages are identified. Notably: 


1. a dense gas stage, leading the emitted gas to move down, eventually touching the ground at a given concentration; this is 
identified as HMP regime 

2. alaminar stage, wherein the fluid moves for a certain distance according to a pure laminar regime 

a dense gas stage consisting of the gas travel at ground level, referred to as BMQ regime 

4. the final neutral (Gaussian) dispersion stage 


po 
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The first three stages have been approached, according to the sequential adoption of validated models and referenced equations, as 
detailed in the following sections. 


2 Results and Discussion 
2.1 HMP Regime 


This regime is described according to the dense gas model presented by Hoot, Meroney and Peterka (1973) and recommended by 
CCPS (1996). It applies to upward-pointing releases but provides results also for horizontal releases. Notably, regarding Figure 2, the 
following resolutive equations are provided. The model has been adopted to assess the outcome of important dense gas dispersion 
scenarios, such as the Jack Rabbit II trial 8 chlorine field experiment (Hanna et al., 2021). 

The maximum initial rise can be calculated according to the following relationship: 


Ah _ 1.32- (“2)F : (2. [ Wo"'Po i (5) 


2:Ro Pa 2-Ro'g'(Po-Pa) 


where: AH is the maximum initial rise, m, R, is the stack/hole radius, m, wo is the gas velocity at the outlet, m/s, u is wind velocity, 
m/s, po is the gas density at the outlet conditions, kg/m’, p4 is the ambient density, kg/m%, g is the acceleration of gravity, m/s”. 


The touch-down distance x, at which the centreline of the dense plume strikes the ground, is given by the analytical relation: 


Xtd _ Wo'U'Po f Ah \3 5 | hs : _ | ? u?-pa y 
2:Ro ~ 2-Ro'9'(Po-Pa) T056 (4) (2 t ud L 2:Ro'g'(Po-Pa) ©) 


where all terms are known, except hs, that is the hole or stack height above ground. 


Finally, the ratio of the maximum concentration at x to the initial concentration C, is: 


— = 2.43. (“) p pae 


o u 2:Ro 
2.2 Laminar Regime 


After impinging the ground, the cloud is assumed to be driven by its true and/or apparent (i.e. cryogenic) density, travelling along 
grade for a given distance, so a dense gas model must be applied. Unfortunately, the most credited model, i.e. Britter McQuaid 
method, is applicable for a concentration range, which falls within the 0.1 and 0.002 fractions of a virtual 100% CO2 initial 
concentration (virtual means that, even if it is the mixture with air, just the CO2 part of the flow is considered). Unfortunately, from 
Xia (100% of the virtual CO2 molar concentration) to 10% of CO», validated simplified models are not available, so the concentration 
profile along this segment is unknown. However, Benintendi (2011), based on the laminar flow of jet flames has adapted the original 
model developed by Aalburg et al. (2005) to the case of laminar flow. This equation can be written as: 


3:CeRgSc 
= FOR 
A EEE (8) 


32er) 


Where xz is the generic distance from the Xm, c:(xı) is the fraction of c(x;a) at x; the value of which falls between 1 (touch down 
distance) and 0.1 (beginning of BMQ regime), Cp is a coefficient is an empirical coefficient ~1, Rg is the Reynolds number, Se is 
Schmidt number. It is worth noting that, in equation 8, all terms are substantially constant, so it describes a hyperbolic function, 
which can be simplified as: 


x, =— (9) 


cL (x) 


where K is therefore constant and must be identified. 
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This function has defined boundary limits, i.e.: cz(0)=1, cx(xz-smo)=0.1, where xz.2muq is the distance from xa where the BMQ regime 
starts. This distance can be calculated by applying the Britter McQuaid method, then, backwards, K is defined in equation 9, so that 
the concentration profile along the laminar regime is fully defined. 


2.3 Britter McQuaid Regime dá 


As stated, the BMQ model (1988) identifies concentrations of heavy X 4/2 
gas on the ground against distances from the release, starting from 0.1 
fraction of the initial concentration up to 0.002 fraction of the same 
initial concentration. It has also been highlighted that the definition of 


the distance from the release at 0.1 allows one to obtain backwards the 
profile of the laminar regime. Figure 3. shows the nomogram 
provided by the authors of the model for continuous releases, with the 


Te- 


i =~ 
ULL-SCALE 


DATA REGION 
definition of the following terms: | 
2 '0=2(PoPa)/pa, m/s”, being po and pa initial plume density and air A 
density respectively, Uref is wind velocity, m/s, qo is the initial 101 č l L v 


E 
u 


ref 


volumetric flow, m?/s, cm is the concentration at distance x, c, is the 
initial concentration 

. N p . Fig. 3 Britter McQuaid nomogram for continuous release. 
Megaris Platform has implemented the mathematical functions 
describing the curves included in the nomogram, both for instantaneous and continuous releases. Once the gas has attained the 0.002 
fraction of the initial concentration, the mixture may likely be 


considered to evolve according to passive dispersion and the Table 4 Inputs and outputs for case study 1. 


cloud will not present any toxic hazard, as the maximum 


HMP Dispersion Model symbol Unit Value 
possible concentration would be 2000 ppm for a pure CO2 Stack height h, = 8 
flow. 
Outlet Temperature To K 223.16 
2.4 Case Studies Outlet radius Ro m 0.3 
Two case studies have been analysed, the former being an Outlet velocity Wo m/s 15 
upward release, and the latter a horizontal release. Wind speed u m/s 1.5 
Ambient Temperature: Ta K 288.16 
2.4.1 Case Study 1-Upward Release Molecular Weight nee shiol eT 
: ; ara ‘ Outlet trati Co / 1 
Including this concentration in equation 4, results from Pr<0, aaa ik aa AA 
so it is not necessary to proceed further, as no toxic conditions Release direction 7 : Vertical 
exist, based on the probability of death obtainable from Table Flow rate Qw kg/s 4.241 
2, which is equal to zero. Air Density Pa kg/m? 1.226 
However, the touchdown concentration can imply some Outlet density ra kg/m? 2.403 
occupational issues. This circumstance should drive the correct f 
Stns ; Jet rise Ah m 9.1 
localisation of plant sections, where peoples are expected to 
stand for long periods of time. Touchdown distanc Su a 41.2 
The values of the touchdown concentration and distance are _ Touchdown concentration Cua ppm 15334 


consistent with those obtained through the most reputable commercial software. 


2.4.2 Case Study 2 - Horizontal Release 
The horizontal release has been resolved, showing a significant concentration Cıa, 62720, at the touchdown distance Xia, 77.2 m, from 


the release point. This value is greater than the IDLH, which is 40000 ppm, according to NIOSH. Consequently, a further analysis of 
the following regimes is required. Input data and findings are shown on Table 5 and in Figure 4. 
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Table 5 Inputs and Outputs for Case Study 2. 
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HMP - Laminar — BMQ Dispersion Model 
HMP Regime 


Item symbol Unit Value Remarks 

Stack height hs m 8 

Outlet Temperature: To K 252.72 

Outlet radius Ro m 0.05 

Outlet velocity Wo m/s 200 

Wind speed m/s 1.5 

Ambient Temperature Ta K 288.16 

Molecular Weight MW g/mol 44.01 

Starting concentration Ca v/v 1 

Release direction - - horizontal 

Mass Flow rate Qw kg/s 1.571 

Air Density Pa kg/m3 1.226 

Outlet density Po kg/m3 2.403 

Jet rise AH m 0.01 

Touchdown distance Xia m TTZ 

Touchdown concentration Cia ppm 62720 

Laminar Regime 

Length of Laminar Regime XL-BMQ m 49.4 Calculated from BMQ 
Constant of Aalburg Equation K m 6272 

Temperature To K 285.71 Heat balance (CO, + Air) 
Duration - - - Continuous Plume 

CO, mass flow rate G kg/s 1.571 

Windspeed at 10 m u m/s 1.5 

Ambient Temperature: Ta K 288.16 

Concentration at the laminar regime start Cu ppm 62720 From HMP 

Concentration at a distance from TD m 0 62720 ppm HMP 
Concentration at a distance from TD m 9.88 31360 ppm Laminar regime 
Concentration at a distance from TD m 19.76 15680 ppm Laminar regime 
Concentration at a distance from TD do 29.64 10453 ppm Laminar regime 
Concentration at a distance from TD m 39.52 7840 ppm Laminar regime 
Concentration at a distance from TD u 49.4 6272 ppm Laminar regime 
BMQ Regime 

Concentration at a distance from TD m 49.4 6272 ppm Britter McQuaid 
Concentration at a distance from TD m 50 6132 ppm Britter McQuaid 
Concentration at a distance from TD m 75 2378 ppm Britter McQuaid 
Concentration at a distance from TD m 100 1257 ppm Britter McQuaid 
Concentration at a distance from TD m 150 652 ppm Britter McQuaid 
Concentration at a distance from TD m 300 211 ppm Britter McQuaid 
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The findings of the two case studies 
confirm the adequacy of the methodology 
to represent the behaviour of a dense gas 
along its evolution from release to the 
start of the Gaussian dispersion. In the 
scientific literature, the various regimes 
are always presented with stand-alone 
methods, so it is not immediate to predict 
the outcome of a release, and, in turn, to 
identify the toxic effects of the resulting 
cloud. The HMP method is shown to work 
very well both for vertical and horizontal 
release. Moreover, the limits of the BMQ, 
which is valid only for concentrations 
starting from the 0.1 fraction threshold 
onward, have been overcome through the 
introduction of a laminar regime equation, 
which allows one to follow the cloud 
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50 100 150 200 250 300 350 


Fig. 4 Profile of CO, concentration for Case Study 2. 


within the more critical concentration bad, as shown in Figure 4. The approach is deemed to represent an innovative method for risk 


assessment of dense gases. 


Conclusions 


CCS is a fundamental process to reduce the impact of carbon dioxide on climate change. The article has presented a very simple but 
comprehensive method, which can support the process safety engineers in CCS studies. The method encompasses more regimes, 
which is infrequent in the normal analysis of the releases and has introduced a new unpublished calculation segment, named Laminar 
Regime, to cover the gap of the concentration profile before the start of the BMQ regime. The calculations have been carried out by 
means of Megaris Platform, a powerful system including data from up to ten thousand chemicals and dozens of calculation models. 

The methodology presented is a valuable tool, even because has the potential to develop the sensitivity capabilities of safety 
engineers, increasing their attitude to govern the analysis scenarios. 


Nomenclature 


aandb  =unitless coefficients 

BMQ =Britter McQuaid 

Cc =the substance concentration 

Ce =initial concentration 

Cx =the initial concentration 

CCPS =Centre for Chemical Process Safety 

C; =a coefficient is an empirical coefficient ~1, 
Cin =the concentration at distance x 

Ah =plume rise 

AH =the maximum initial rise 

DTL =Dangerous Toxic Load 

g =the acceleration of gravity 

HMP =Hoot Meroney Peterka 

hs, =stack or hole height 

HSE =Health and Safety Executive 

IDLH =Immediately Dangerous to Life and Health 
JT =Joule Thomson 

K =constant of Aarburg equation 

LOC =Loss Of Containment 

NIOSH =National Institute for Occupational Safety and Health 
do =the initial volumetric flow 

QRA =Quantitative Risk Assessment 

Ro -the stack/hole radius 

Rg =Reynolds number 

Pa =the gas density at the outlet conditions 

Pa =the ambient density 

Se =Schmidt number 


SLOD =Significant Likelihood of Death 


64 


Jordanian Journal of Engineering and Chemical Industries (JJEC! 


SLOT =Specified Level of Toxicity 
STEL =Short-Term Exposure Limit 
To. =initial Tempertaure of CO, flow 
t =the exposure time 

TD =TouchDown 

TLV =Threshold Limit Value 
TWA =Time-Weighted Average 

u =wind velocity 

Uref =wind velocity 

Wcoz =C0, Mass Flow 

Wo =the gas velocity at the outlet 
Xid =touch down distance 

XL =generic distance from Xia 
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